Dynamics of a quantum phase transition with decoherence: 
the quantum Ising chain in a static spin environment 
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We consider a linear quench from the paramagnetic to ferromagnetic phase in the quantum Ising 
chain interacting with a static spin environment. Both decoherence from the environment and non- 
adiabaticity of the evolution near a critical point excite the system from the final ferromagnetic 
ground state. For weak decoherence and relatively fast quenches the excitation energy, proportional 
to the number of kinks in the final state, decays like an inverse square root of a quench time, 
but slow transitions or strong decoherence make it decay in a much slower logarithmic way. We 
also find that fidelity between the final ferromagnetic ground state and a final state after a quench 
decays exponentially with a size of a chain, with a decay rate proportional to average density 
of excited kinks, and a proportionality factor evolving from 1.3 for weak decoherence and fast 
quenches to approximately I for slow transitions or strong decoherence. Simultaneously, correlations 
between kinks randomly distributed along the chain evolve from a near-crystalline anti-bunching to 
a Poissonian distribution of kinks in a number of isolated Anderson localization centers randomly 
scattered along the chain. 

PACS numbers: 75.10.Pq, 03.65.-w, 64.60.-i, 73.43.Nq 
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I. INTRODUCTION 

Phase transition is a fundamental change in the state 
of a system when one of its parameters passes through 
the critical point. In a second order phase transition, the 
fundamental change is continuous and the critical point is 
characterized by divergences in the correlation length and 
in the relaxation time. This critical slowing down implies 
that no matter how slowly a system is driven through the 
transition, its evolution cannot be adiabatic close to the 
critical point. As a result, ordering of the state after the 
transition is not perfect: it is a mosaic of ordered domains 
whose finite size £ depends on the rate of the transition. 
This scenario was first described in the cosmological con- 
text by Kibble [J who appealed to relativistic casuality 
to set the size of the domains. The dynamical mechanism 
relevant for second order phase transitions was proposed 
by Zurek 0. It is based on the universality of critical 
slowing down, and leads to a prediction that average size 
£ of the ordered domains scales with the transition time 
tq as 



£ 



(1) 



where v and z are critical exponents. The Kibble- 
Zurek mechanism (KZM) for second order thermody- 
namic phase transitions was confirmed by numerical sim- 
ulations of the time-dependent Ginzburg-Landau model 
[H and successfully tested by experiments in liquid crys- 
tals i_0| , superfluid helium 3 [5[ , both high-T c [y] and low- 
T c [7[ superconductors, and even in non-equilibrium sys- 
tems With the exception of superfluid 4 He - where 
the early detection of copious defect formation was 
subsequently attributed to vorticity inadvertently intro- 
duced by stirring 10], and the situation remains un- 
clear - experimental results are consistent with KZM, al- 



though more experimental work is clearly needed to allow 
for more stringent experimental tests of KZM. Quite re- 
cently, a new experiment was reported in Ref. where 
they observe, for the first time, spontaneous appearance 
of vorticity during Bose-Einstein condensation driven by 
evaporative cooling confirming KZM predictions in Ref. 

m. 

The Kibblc-Zurek mechanism is thus a universal the- 
ory of the dynamics of second order phase transitions 
whose applications range from the low temperature Bose- 
Einstein condensation (BEC) to the ultra high temper- 
ature transitions in the grand unified theories of high 
energy physics. However, the zero temperature quan- 
tum limit remained unexplored until very recently and 
quantum phase transitions are in many respects quali- 
tatively different from transitions at finite temperature. 
Most importantly time evolution is unitary, so there is 
no damping, and there are no thermal fluctuations that 
initiate symmetry breaking in KZM. The recent progress 
on dynamical quantum phase transitions is mostly the- 
oretical^*; Refs. Ill El EE EE EE El, El, Eft HH, 
HI, HI Hi, HE HE and, for an example of a disor- 
dered quantum system, Ref. [28[, but there is already 
one significant exception: the experiment in Ref. [291 ] on 
a transition from paramagnetic to ferromagnetic phase 
in a dipolar BEC. Generic outcome of that experiment is 
a mosaic of finite-size ferromagnetic domains whose ori- 
gin was attributed to the Kibble-Zurek mechanism. This 
explanation is further supported by theory in Ref. [2lj |. 

A majority of theoretical work was devoted to the pro- 
totypical exactly solvable quantum Ising chain [IE El, 
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driven by a linear quench 

git) 
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from g — oo to g — i.e. across the quantum phase tran- 
sition from paramagnet to ferromagnet at g c = 1. Since 
the system is gapless at g c , when tq > 1 the evolution 
becomes non-adiabatic at 



(4) 



see Ref. [15]. The freezeout at g is the closer to g c the 
slower is the transition. The Ising critical exponents v = 
1 and z = 1 determine an average size of ferromagnetic 
domains at g = proportional to the correlation length 
at g: 



t „ 1/2 

or final density of kinks (domain walls) at g = 



r 1 



-1/2 



(5) 



(6) 



which is proportional to excitation energy density. Other 
properties, like spin-spin correlation functions flS. l20l| or 
entropy of entanglement between a block of consecutive 
spins and the rest of the chain [2C| were analyzed in some 
detail and they all turned out to be determined by £. 

While the quench in the isolated system ([2]) seems to 
be well analyzed, relatively little is still known about dy- 
namical transitions in open quantum systems subject to 
interaction with environment. Significant progress was 
made in Ref. [30] in a "global" case when interaction 
between the system S with the Hamiltonian ^ and its 
environment E is described by V = R (J2 n °n)' wnere 
is a hermitian operator of the environment. This global 
model is solvable thanks to its translational invariance 
and its solutions indicate that decoherence is increasing 
density of excited kinks as compared to an isolated sys- 
tem. A local model, with the system ([^coupled to an 
Ohmic heat bath, was analyzed in Ref. [311 ] distinguish- 
ing between different regimes of parameters where defect 
production is dominated either by KZM or external heat- 
ing. In this paper, we propose a quite realistic, but still 
solvable, model of local zero temperature decoherence 
from a static environment. Its solution predicts dramatic 
increase in the density d of excited kinks as compared to 
an isolated system with d decaying as only a logarithmic 
function of tq. 

Motivation for this study is twofold. It comes both 
from condensed matter physics, where it is virtually im- 
possible to isolate a system from its environment, and 
adiabatic quantum computation, where a system is ini- 
tially prepared in a simple ground state of a simple initial 
Hamiltonian Hq and then it is driven adiabatically to a 
final Hamiltonian Hi whose non-trivial ground state is 
the desired solution of a complex computational prob- 
lem. The computation is complicated by a quantum crit- 
ical point somewhere on the way from Hq to H\ which 



can make the adiabaticity problematic, but see Refs. 
[2^. [27| for methods how to circumvent this problem. The 
Ising chain ^ is a toy model of adiabatic quantum com- 
puter with the final (trivial) ferromagnetic ground state 
at g — playing the role of the desired "non-trivial" 
ground state. When the "computer" is isolated from en- 
vironment, then Eq. ^ implies that the minimal "com- 
putation time" tq required to keep the evolution adia- 
batic or, equivalently, to make f ^> N is 



isolated 



N 



(7) 



The "isolated" computation problem is polynomial in N. 
In contrast, in our model of decoherence similar argument 
predicts 



open r~J gVJV 



(8) 



which is non-polynomial in N. 



II. ISING CHAIN IN STATIC SPIN BATH 

In this paper we couple the Ising chain ^ to an envi- 
ronment E of M spins through the interaction 



N M 



(9) 



Here r m 's are Pauli matrices of environmental spins. The 
spins are static, with He = 0, and the total Hamiltonian 
is just H = H s + V. 

Initially at t — > —00 the system is in the ground state 
|0 g _»oo) of the pure Ising chain with all spins polar- 
ized along x. This assumption is self-consistent in our 
open system because large initial energy gap of 2g makes 
the influence of the static environment so negligible that 
the initial states of S and E can be assumed uncorre- 
cted: ps+E — Ps ® Pe with p s = |0 g ^ oo )(0 g ^ oo | and 
the environment is initially in a pure state 



E 



si 



S M 



(10) 



,S A f = -l, + l 



Here T? n \s m ) = s m \s m/ 



After evolution for time At reduced density matrix of 
the system ps = Tr ePs+e becomes 



Ps(At) = ^|c s f U(At,s) lOooXOool U\At,s) 



= U{At 7 s) lOooXOool W(At,S) 
Here s = (s x , s M ), 



U{At, s) = 7 exp 



i / dt' H(t',s) 
Jo 



(11) 



(12) 
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and 



JY 



H(t, s) = -J2{ [<?(*) + r„] < + « +1 } (is) 



with random magnetic fields of the static environment 



M 



(14) 



The overline in Eq. (jlip is an average over s with prob- 
ability distribution \cg\ 2 , but it can also be interpreted 
as an average over random "disorder" field T n . pg(At) is 
an average over states U (At, s) |0oo) obtained in quenches 
with different disordered Hamiltonians (fTB")) . In this way, 
our original problem of a quench in the open pure Ising 
model ([2]) is mapped to an average over quenches in the 
isolated random Ising model (|13| . 

In the following, rather than struggle with the problem 
in its full generality, we assume that each spin of the 
environment couples to only one spin of the system or, 
in other words, each spin of the system has its own local 
environment. Consequently, T m and T n are statistically 
independent when m ^ n. We also assume that each T n 
has the same Gaussian probability distribution 



/(r) 



V2~t 



(15) 



where a is strength of disorder/decoherence. 

The Hamiltonian (fTB")) belongs to the universality class 
of the well known random Ising chain [32, HH . It has a 
continuous quantum critical point at g c when 



dT f(T) ln|.g c + r| = . (16) 



g c depends on a as shown in Fig. [T|A There is no critical 
point for a > 1.887 when the disorder is too strong. No 
matter how weak a is, renormalization group transfor- 
mations drive the model (|13p towards an infinite disorder 
fixed point with different critical exponents than in the 
pure Ising chain: the random chain has v = 2 instead of 
v = 1 and z — > oo instead of z = 1 .'!•'! . A straightfor- 
ward application of the standard KZ formula |T]) gives 
£ ~ 1 i.e. a domain size that does not depend on the 
quench time tq. A more careful argument in Ref. [28| . 
going back to the basics of KZM, predicts a logarithmic 
dependence 



In [ln(ar Q )] 



(17) 



with a non-universal a ~ 1. This equation leads to the 
estimate in Eq. ©. Since Eq. (fT7]) is based on the 
universality class alone, it is valid for any model of this 
class when tq is long enough for the quench to become 
non-adiabatic close enough to g c to be affected by disor- 
der. The estimate (JTTJ) was confirmed by numerics in the 



model of Ref. [28j , where it was also found that for weak 
disorder and relatively fast tq one recovers Eqs. (|5I6[) as 
in the pure model Eq. ©. 

In the present effective model (fTB")) , it is simple to esti- 
mate how slow a quench needs to be for Eqs. (|5l6p and, 
more importantly, Eq. to be not valid. Assuming 
that influence of T n is negligible, evolution becomes non- 
adiabatic at a field g in Eq. This assumption is not 
self-consistent when the remaining distance from g to g Cl 

A —1/2 

9 ~ 9c — t q , is less than the strength a of disorder 
field r n , or equivalcntly 



> 1 



(18) 



Thus, no matter how weak the decoherence is, its influ- 
ence is not negligible when the transition is slow enough: 
tq 3> er~ 2 . In consequence, there is a maximal number 
of qubits 



N < 



1 



(19) 



which can be simulated with polynomial efficiency, com- 
pare Eq. © and ([15]). 

In the next Section, we review static properties of the 
random Ising model (|13p many of which, we believe, are 
described in this form for the first time. 



III. RANDOM ISING CHAIN 

Here we assume for convenience that N is even and fol- 
lowing Refs. [28|, HH make the Jordan- Wigner transfor- 
mation a x n = \- 2clc n and a z n = - (c„ + 4) ]J m<n (l - 
2c\ n c m ) introducing spinless fermionic operators c„. The 
Hamiltonian M becomes H = P+H+P+ + P-H~P~ 



where P 



2cJ,c n ) 



on subspaces with even (+) and odd 
quasiparticles and 



are projectors 
) numbers of o 



JV 



71=1 

h.c. 



are quadratic Hamiltonians. Here 
9n = g + r„ 



(20) 



(21) 



for short. The c„'s in H~ satisfy periodic boundary 
conditions c^v+i = ci, but the c„'s in H + must obey 
c/v + i = — c\, what we call "antiperiodic" boundary con- 
ditions. 

The parity of the number of c-quasiparticles is a good 
quantum number and the ground state has even par- 
ity for any value of g. Assuming that the quench be- 
gins in the ground state we can confine to the sub- 
space of even parity. In this subspace the quadratic 
H + is diagonalized by a Bogoliubov transformation c„ = 



4 




1 

I 1 E 


_ C r at g=0, a=0.8 - 
— P at g=0, a=0.8 " 







FIG. 1: In A the critical g c in Eq. (|16p is shown as a function 
of a in Eq. (115[) . In B density of kinks in a ground state of the 
random chain Eq. (|13[1 is shown as a function of g for N = 512 
and different cr's. In C we show correlation coefficient c in 
Eq. (|31|l as a function of g and a. When a < 1, c f» 0.5 
in the ferromagnetic phase below c/ c ~ 1 and c > 0.5 in the 
paramgnetic phase above g c . In D we show correlator C r in 
Eq. (|36[1 between two kinks in a Cooper pair in the pure Ising 
chain (cr = 0). C r is localized below <? c = 1 and delocalized 
above. In E both CV and P n in Eq. (|37[) are shown in the 
ferromagnetic phase at g = and a = 0.8. Here iV = 512 in 
both D and E. 



Z^m=l \ u nm"/r) 



7jJ. The index m numerates (Bo- 



goliubov) eigenmodes of the stationary Bogoliubov-de 
Gennes equations 



dui 



ill 



2m? 



with uj m > 0. Here we define u^ m = 
sume the anti-periodic boundary conditions: 



(22) 



± v„m and as- 



l N+l,m ~ 

l N,m- Tne eigenstates (u nm ,v nm ), 
normalized so that En (|w„ m | 2 + |iw| 2 ) = 1, define 
quasiparticle operators 7 m = u* m c„ + v *cl. After 
the Bogoliubov transformation the Hamiltonian becomes 



H+ = E 



AT 



idm'Ym ~ I) wmc h is a simple- looking 



sum of quasiparticles. However, thanks to the projection 
P+ H + P+ only states with even numbers of quasipar- 
ticles belong to the spectrum of H . 

When g = a = the final system Hamiltonian ([2|) has 
N degenerate quasiparticles with oj = 2. We choose an 
orthonormal basis 



,(0) 



,(0) 



(Sn+l 
(^n+l,m 



(23) 
(24) 



antiperiodic in n. These eigenmodes define Bogoliubov 
quasiparticles 7^? which are simply kinks localized at 
bonds between sites m and m+1. The kinks are related to 



quasiparticles {u nm ,v nm ) at finite {g,T) by a Bogoliubov 
transformation 



n, - TP J® 4- V* -v (0)t 
7a — u ba 7 fc + v ba 7 & , 



where 



u ba = J2 

n 



(0) 



(0), 



(0), 



(25) 

(26) 
(27) 



leading, for example, to a simple expression for average 
kink density in the ground state |0 S: r) of the Hamiltonian 
(HI 



d{g 7 a) 



^(0 s ,rlE^ 0)t 7^|0,r) 



TrVW 
N 



, (28) 



shown in Fig. [TJ3. It is finite at any g > 0, both in the fer- 
romagnetic (g < g c ) and paramagnetic phase (<? > g c ), 
but, as we will see, kink-kink correlations are qualita- 
tively different in the two phases. 

These correlations can be indirectly probed by average 
fidelity F between the final ground state |0) of Hs at 
.9 = (state with no kinks 7^) and ground states \g, T n ) 
at finite g or cr: 



10 



(29) 



Here Z — V*(U*) 1 and M is a normalization factor. 
The average fidelity is 



F 



<0|0 9 ,r)(0 ff , r |0) 



(N/2 
n=0 v ' 



Z n \0) 



Det(l + ^^) 

„(°)t„,(o)t 



-1/2 



Here Z = § E Qb Z abl ^ lb ^ n for short. 

We found that the fidelity is exponential in N 



F(g,a,N) ~ (1 - c d) 



N 



(30) 



(31) 



when F <SC 1. Here g?(<?, cr) is the density of kinks in Eq. 
(|28|) and Fig. QJ3. The coefficient c(g, a) is shown in Fig. 
[Tp. For weak disorder, when cr <C 1, we have c ss ^ when 
g < (?c with c increasing towards 1 when g c/ c or a 3> 1 . 
These two limits can be explained as follows. 

When the magnetic fields g n = T n + g in the effective 
Hamiltonian (|13p are strong, because either g or cr or 
both are strong, then in any ground state |0 gj r) all spins 
are polarized along cr^. Fidelity to the a z -ferromagnet 
is F = 1/2 N and density of kinks is d = 1/2, hence 
F = (1 — cd) N with c — 1. This c is consistent with the 
data shown in Fig. [Tp when g or a are strong. 



5 



In the opposite limit of weak magnetic fields \g„ 
\g + T n \ <C 1, the ground states are 



N 

n 

n=l 



(32) 



1 + % 
x T 16 



Their fidelity to the cr z -ferromagnet is 



F 



N 

n 



= 1 1 + 16 



N 



(33) 



when F<1. Here d = <C 1 is small density of kinks 
and c = 1/2. This c is consistent with the data in Fig. 
[Tp when cr<l and g < g c ~ 1. 

It is interesting to interpret the widely different values 
of i < c < 1 in terms of a simple Poissonian model where 
each of N bonds is either excited (with probability d cxc ) 
or not excited (with probability 1 — <i e xc) independently 
of other bonds. Here d cxc is average density of excita- 
tions. The fidelity is a probability that none of the N 
independent bonds is excited 



F = ( 1 - de 



(34) 



Comparing Eqs. (|31l34[) we obtain density of independent 
excitations 



c d 



(35) 



We can conclude that c = d cyLC /d measures correlations 
between kinks: c < 1 means bunching and an eventual 
c > 1 would mean anti-bunching of kinks randomly dis- 
tributed along the spin chain. 

This simple interpretation of c follows from the fact 
that any ground state |0 ffi r) is a BCS state of kinks 7W, 
see Eq. (|29|) where Z a b is a wave-function for a (Cooper) 
pair of kinks. Depending on g and cr, this BCS state 
can be either a condensate of tightly bound Cooper pairs 
of kinks (with c = 1/2), or just a weakly coupled BCS 
state with pairing (anti-) correlations manifesting as the 
(anti-) bunching of kinks. 

The BCS state is particularly simple at weak magnetic 
fields g n = g + T n when, crudely speaking, the ground 
state |0 Si r) m Eq. (f3"2"|) is approximately the ferromagnet 
I T1T2 ••• In) but with occasional spins reversed to | J, n ) 
by weak magnetic fields g n ^n- Each reversed spin is a 
tightly bound Cooper pair of two kinks sitting on nearest 
neighbor bonds with a relative distance of r = 1 lattice 
sites between the two kinks. No wonder that d eX c = d/2 
in this case. 

As we could see in Fig. [T]C, the picture with tightly 
bound Cooper pairs is accurate for weak disorder in the 
ferromagnetic phase, but with g increasing above g c the 
Cooper pairs begin to dissociate into free kinks and an- 
tikinks and c increases well above \. This dissociation 
can be clearly seen in Fig. [Tp where we show a correla- 
tion function C r 



\z. 



m-\-r,m \ 



(36) 



between two kinks making a Cooper pair. Here Z a ^ is the 
unnormalized wavefunction of a Cooper pair of two kinks 
in the BCS state (|29|) and r is a relative distance between 
the two kinks. The C r shown in Fig. QJ) is localized 
when g <= g c (tight Cooper pairs) and delocalized when 
g > g c (dissociated pairs). 

A similar C' r is shown in Fig. [TJi! for a finite a — 0.8 in 
the ferromagnetic phase at g = 0. In addition, the same 
figure shows reduced probability distribution 



Pa 



\z n 



(37) 



for a single kink. Unlike in the pure case, this P„ is not 
uniform and shows Anderson localization. 

In the ferromagnetic phase below g c all kinks are bound 
into tight Cooper pairs. The small density of Cooper 
pairs, each of them a tightly bound pair of kink and an- 
tikink, does not destroy long range ferromagnetic order 
in a z . With g increasing above g c the transition to the 
paramagnetic phase takes place when the Cooper pairs 
begin to dissociate into free kinks and antikinks which 
scramble long-range ferromagnetic correlations. 

In the context of adiabatic quantum computation there 
are two generally accepted measures of how far the final 
state ps(0) from the desired final ground state |0) is. 
One is energy of excitation of the final state above the 
energy of the desired final ground state, and the other 
is fidelity between the final state and the final ground 
state. Both quantities are tractable in our model, where 
the excitation energy is proportional to the number of 
kinks. We study them in the following two Sections. 



IV. DENSITY OF KINKS AFTER A QUENCH 

In our simulations of a quench the system is initially 
prepared in its ground state at a large initial value of 
g 1 i.e. in a Bogoliubov vacuum state for quasiparti- 
cles at an initial g ^> 1. As the magnetic field is being 
turned off to zero across <? c , the state of the system \tp(t)) 
is getting excited from its adiabatic ground state. How- 
ever, in a similar way as in Refs. [l6|, [28[, we use the 
Heisenberg picture where the state remains a vacuum for 
quasiparticle operators 



In 



(38) 



with the Bogoliubov modes u nm (t) and v nm (t) solving 
time-dependent Bogoliubov-de Gennes equations 



duz 



dt 



1,771 



(39) 



Initially each mode is a positive frequency eigenmode of 
the stationary Eqs. ([22|) . Equations (|39|) were integrated 
by a split-step method for each realization of r„. Their 
solutions u^ m (0) at the final g = determine final states 
IV'(O)) whose average over disorder gives final density ma- 
trix p s (0) = |^(0)) (^(0) I . For each realization of r„, a 



g(t) 

FIG. 2: Density of kinks during a quench as a function of 
g(i). Here each plot is a single realization with a — 0.1 on 
a lattice of N — 512 sites. For the early large g the density 
follows the density d(g,a) in a ground state |0 9 ,r), but as the 
quench is approaching the critical point g c ~ 1 the evolution 
becomes non-adiabatic and the density gets excited above its 
ground state level. The slower is the quench, the closer to the 
critical point the non-adiabatic stage begins and the quench 
is probing the critical point more closely. 
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FIG. 3: In A final density of kinks at g = is shown as a 
function of tq on a N = 512 chain. For large tq the density 
tends to saturate at the average density of kinks d 9 =o,<r in 
the ground state of the random Hamiltonian (|13[) shown in 
Fig. (QJB. In B we show difference Sd — d — d g= o, a which is 
density of kinks excited above the ground state of the random 
Hamiltonian at g = 0. This difference can be attributed to 
the non-adiabaticity of the transition described by KZM. 



final state IV'(O)) is a vacuum for 7 m 's which are related 
to the kinks jm by the transformation Eq. ([25]) . Final 
densities of kinks follow from Eq. (f2"5| . 

For the system Hamiltonian @ at the final g — the 
final excitation energy is simply twice the final number of 
kinks excited in the desired ferromagnetic ground state. 
In Figure [5] we show density of kinks d as a function of 
g(t). For the early large g the density follows the density 
d(g,a) in a ground state |0 Sl r), which is also shown in 
Fig. [T]B, but as the quench is approaching the critical 
point g c ps 1 the evolution becomes non-adiabatic and the 
density d is excited above its ground state level d(g,a). 
The slower is the quench, the closer to the critical point 
the non-adiabatic stage begins and the quench is probing 
the critical point more closely. 

The final density of kinks at g = is shown in Fig. [3)^. 
For large tq, the density tends to saturate at d{g — 0, a) 
shown in Fig. (HJB i.e. the density of kinks in the 
ground state of the random Hamiltonian Q13p . In Fig. 
|3]B we show a difference Sd = d — d(g = 0, a) which 
can be attributed to the non-adiabaticity of the tran- 
sition described by KZM. If Sd = o:Tq, then in the 
log-log plot of Fig. [3J3 we would see a line log 10 Sd = 
log 10 a + w(log 10 tq), but this is not the case when 
tq ^> cr~ 2 as in Eq. (fT5|) . At best we can think of a 
local slope w(tq) which can be estimated by fitting to 
pairs of nearest neighbor data points. In the legend we 
give ranges of local slopes w obtained for different cr's 
(with error bars on their last digits). For weak a and 
small tq the slope w is close to the —1/2 characteristic 
for the pure model fast quenches, when tq <§; <r~ 2 , 



become non-adiabatic far enough from the critical point 
not to see any effect of weak disorder. At stronger a or 
longer tq the local slopes are less steep and for a fixed a 
they become less steep with increasing tq. For example, 
at the strongest a = 0.8 the local slope falls to a mere 
| w | = 0.04 for the longest tq. These observations are 
consistent with the predicted logarithmic dependence of 
the dynamical correlation length £ in Eq. (fTT)) . 



V. FIDELITY AND CORRELATIONS AFTER A 
QUENCH 

In the context of adiabatic quantum computation, it is 
important to know how close the final state ps(0) to the 
desired ground state of the final Hamiltonian Hs is. The 
closeness is measured by fidelity 



\ PS (0)\0) = (0|V(0))(^(0)|0) (40) 



F 



given by Eq. ([30|) . 

Without decoherence, or in a pure Ising chain with 
(7 = 0, the fidelity can be obtained analytically from 
the exact solution in Refs. [l|| [2(|. F is a probability 
that not a single pair of 7'°'-quasiparticles with opposite 
quasimomenta (k, —k) is excited after a quench, F = 
rifc>o( 1 ~Pk)- Here p k ~ exp(-~2TTT Q k 2 ) when t q > 1. 
When Nd >lwe obtain 



kiF 



Nd 



f™ds In 



1 



Jo ds e ~ 



-l.SNd . (41) 
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FIG. 4: In A and C final kink density d for different tq as 
functions of lattice size N for a — 0.1 (panel A) and a — 0.8 
(panel C). In B and D fidelity F for a = 0.1 (panel B) and 
a — 0.8 (panel D). With inreasing N the density saturates 
at an asymptotic d(jQ,a) when N 3> l/d(jQ,a). In the 
same asymptotic regime the fidelity becomes exponential in 
N . Here the averages are taken over Nr realizations with 
N R N > 2048. 
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FIG. 5: Correlation coefficients c for different tq and a ob- 
tained by linear fits to the exponential tails of fidelity in Figs. 
HP and D. 



Here d = f*^ ^Pk — l/27r-y/2rQ. F is exponential in 
N as in the static case Eq. ([31]) . Given that d <C 1 
for tq 1 we obtain that c ~ 1.3 > 1 implying anti- 
bunching of kinks. 

The anti-bunching can also be seen in the state | -0(0)) 
after a quench which is a BCS state of kinks in Eq. (f2U)) . 
Figure [6]A. shows a probability distribution P n for a kink 
in Eq. (|37|) and the correlation function C r between two 
kinks in a Cooper pair in Eq. (|36[) after a quench with 
(7 = and tq = 16. We have Co = because the kinks 
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FIG. 6: In A a correlation function CV between two kinks in 
a Cooper pair and probability distribution P n for a kink on a 
iV = 512 periodic lattice after a quench with a — and tq — 
16. In B, C, D, and E we show typical final CV and P„ after 
quenches with a — 0.1, 0.8 and tq — 16, 1024. Here all C r and 
P n for a given a come from the same realization of disorder 
r n . Both C r and P n show Anderson localization. When tq or 
a are large, then CV is localized around those r's which can be 
identified as distances (modulo periodic boundary conditions) 
between the localization centers in the corresponding P n , with 
the exception of r ~ which is avoided because two fermionic 
kinks do not like to choose the same localization center. 



are fermions, there is a broad maximum in C r in the 
range \r\ = 20... 40, and a flat distribution for |r| > 
40 when the Cooper pair is dissociated. If C r were flat 
everywhere (except r = 0), then c = 1, but the broad 
maximum means that even when the kinks get close to 
each other they prefer to keep a safe distance in the range 
20 . . . 40. This short range repulsion is consistent with the 
anti-bunching observed in Eq. (|41|) where c = 1.3 > 1. 

The same anti-bunching can also be seen in the ferro- 
magnetic correlation function after a quench at g = 



cr - a , 



exp (-l.55.Rd) cos (2.95-Rd - ip) (42) 



accurate when 1 <C R <C ^/tq log tq, see Ref. [2Qj. The 
oscillatory cosine factor means that kinks tend to order 
into a crystal lattice. This very imperfect crystalline or- 
der implies the anti-bunching seen in c = 1.3 > 1. We 
can conclude that in the state after a quench with a = 
there is a similar connection between c, C r , and ferro- 
magnetic correlations as in the static case. 

For quenches with a > we also find an exponential 
tail 



F( TQ ,a)~(l-cd) 



N 



(43) 



when F <C 1, compare panels B and D in Fig. [U Here 
d(TQ,a) is an asymptotic value of average kink density 
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FIG. 7: In B and C convolutions in Eq. $7} corresponding 
to C r 's in Figs. [6)3 and D obtained from distributions P„ in 
Figs. [Bp and E. Except for tq — 16, a — 0.1 (black in panels 
B), the convolutions are very close to their corresponding CV's 
everywhere apart from r close to (and N). 



obtained for a sufficiently large lattice size N, such that 
NcI(tq,<t) ^> 1, see panels A and C in the same Fig. 
[U In contrast, when N <C 1/c£(tq, er), then a finite gap 
at g = g c results in an adiabatic transition and a final 
density of kinks which is less than the asymptotic cI(tq , a) 
for large N. 

Fits to the exponential tails of fidelity in Figs. 0)3 and 
D give correlation coefficients c for different tq and a 
which are shown in Fig. [5) For comparison, we also show 
in the same figure c at a = which saturates at c = 1.3 
when tq 3> 1. The results suggest that for a strong a or 
large tq, when even a weak a has strong effect, the corre- 
laton coefficient decays towards c ~ 1 suggesting random 
Poissonian trains of kinks and exponential ferromagnetic 
correlation functions. 

This picture is corroborated by typical correlation 
functions C r and probability distributions P n shown in 
panels B, C, D, and E of Fig. [5] Here it is clear, espe- 
cially for the larger tq — 1024 (green) or the stronger 
a = 0.8 (panels C and E), that a kink gets localized in 
isolated Anderson localization centers. A close inspec- 
tion of corresponding CV's reveals that C r is localized at 
r's equal to distances between the localization centers in 
P n . There is only one exception from this rule: C r is neg- 
ligible when r fts or N because, apparetly, a (Cooper) 
pair of fermionic kinks avoids being trapped in the same 
localization center. It seems that each kink in a Cooper 
pair chooses its localization center at random, indepen- 
dently of the other kink, except for avoiding the same 
localization center. A quantitative proof of this simple 
picture is provided in Figs. 03 and C where we plot a 
convolution 



PPr 



N 

E 

n=l 



P P 



(44) 



If the two kinks in a Cooper pair were independent, 
then this convolution would reconstruct the correspond- 
ing CV's in Figs. [BJ3 and C. Comparing panels B and 
C in the two figures we see that it does, with the ex- 
pected exceptions when r « 0, N and in the case of a 
fast quench (tq = 16) at weak <r = 0.1 (black in Fig. 
HJ3). We can conlude that, with some idealization, for 
long tq or strong cr, final kinks are distributed as if each 
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FIG. 8: In A CV and P„ in a ground state of the random 
Hamiltonian (|13[) at a — 0.8 and g = 1.3 just above g c . Here 
P n shows similar fragmentation into isolated Anderson local- 
ization centers as in the final states after a quench with long 
tq or large a in Figs. [5] In B we show both the CV from 
panel A and a convolution PP„ in Eq. (|44[) obtained form 
the P n in panel A. The two plots are identical, except for 
r = 0, demonstrating independence of kinks. 



kink were choosing at random one of the isolated Ander- 
son localization centers randomly distributed along the 
chain. This approximate Poissonian model is consistent 
with the observed c ~ 1 . 

When tqo 1 1, then the final states at g = are 
qualitatively different from the ground state of the ran- 
dom Hamiltonian (p~3|) at g = 0: the (green) fragmented 
C r 's in Figs. [SJ3 and C are qualitatively different from 
the (black) C r in Fig. [TJH localized around r = 0. In 
other words, all kinks in the ground state, of density 
d(g = 0,er), are tightly bound into Cooper pairs which 
do not destroy ferromagnetic long range order, while all 
kinks in a final state, of density d, contribute to exponen- 
tially decaying ferromagnetic correlations even though 
for slow quenches the difference 5d = d — d(0, a), whose 
origin can be attributed to non-adiabaticity, is small as 
compared to d(0, cr). 

These striking properties of the final states after a 
quench are inherited from the ground state of the ran- 
dom Hamiltonian (| 13[) just above g c . Its properties are 
relevant here because, as we could see in Fig. [21 in a 
slow quench a state \tp(t)) follows an adiabatic ground 
state of the random Hamiltonian until a g just above g c 
when, in accordance with KZM, the evolution becomes 
non-adiabatic. In Fig. [51 we show some properties of the 
ground state just above g c where P n is fragmented into 
isolated localization centers, and C r is identical with PP r 
everywhere except r = 0. Apparently, the following evo- 
lution from g to g — does not change this qualitative 
picture. 



Conclusion 

A static spin environment increases non-adiabaticity of 
the transition in a dramatic way: density of quasiparti- 
cles (kinks) decays no longer as a power of the transition 
time but in a much slower logarithmic way. This means, 
in the context of adiabatic quantum computation, that 
coupling to a static environment may transform a polyno- 
mial computational problem into a non-polynomial one. 
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Fidelity between a final state after a quench and the 
desired final ground state decays exponentially with a 
chain size. The rate of this decay is equal to the density 
of kinks times a correlation coefficient equal, for fast tran- 
sitions and weak decoherence, to 1.3 and, for slow tran- 
sitions or strong decoherence, close to 1. Corresponding 
kink-kink correlations are, respectively, anti-bunching in 
a near-crystalline ordering and a simple Poissonian distri- 
bution of kinks in isolated Anderson localization centers 
randomly distributed along a chain. 
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